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We use a hydrodynamics theory to discuss the response of a viscous flux-line liquid 
to an ac perturbation applied at the surface of the sample. The theory incorporates 
viscoelastic effects and describes the crossover between liquid-like and solid-like response 
of the vortex array as the frequency of the perturbation increases. A large viscosity from 
flux-line interactions and entanglement leads to viscous screening of surface fields. As a 
result, two frequency-dependent length scales are needed to describe the penetration of 
an ac field. For large viscosities the imaginary part of the ac permeability can exihibit, 
in addition to the well-know peak associated with flux diffusion across the sample, a new 
low-frequency peak corresponding to the transition from solid-like to liquid-like behavior. 



3/94 



PACS numbers: 74.60.Ec, 74.60.Ge, 74.40.+k. 



1. Introduction 

The study of the response of high-temperature (HTC) superconductors in the mixed 
state to alternating magnetic fields and transport currents provides direct information on 
the dynamics of magnetic flux lines in these materials Such ac measurements 

can be performed by superimposing a small alternating magnetic field to the constant 
applied field, H. Alternatively, one can measure the change in resonance frequency of the 
superconducting crystal undergoing mechanical vibrations ]5|,[| . Due to the pinning of the 
flux lines to the crystal lattice, the periodic tilting of the superconductor is equivalent to 
the application of an ac field normal to H. 

AC techniques are also used to determine the irreversibility line, Ti rr (H), in disor- 
dered HTC superconductors J7|. In magnetization measurements this is defined as the 
locus in the (H : T) phase diagram where the magnetization for field-cooled samples dif- 
fers from the result obtained in samples cooled initially at zero field. This criterion is 
experimentally ambiguous since it depends on the rate at which data are taken. For this 
reason ac techniques, such as ac permeability, ac transport measurements, ac field pene- 
tration, mechanical oscillator |5|,[| , and others are often preferable. When cooling through 
the irreversibility line during an ac experiment the dissipation goes through a maximum. 
Experimentally the location of this maximum is very close to Ti rr (H) and is used as an 
alternative definition of Ti rr . The location of Ti rr (H) has great technological as well as 
intellectual interset since it generally coincides with the line below which the resistivity 
becomes unmeasurably small. 

Much theoretical work has been done to investigate the high-frequency electro dynamic 
response of a vortex solid pinned by weak impurity disorder |8|,|lT] • When the perturbing 
field is step-like in time, the perturbation initially interacts with the flux line only at the 
surface of the sample. The field, transport current or mechanical tilt generate a surface 
current that flows in a layer of thickness of the order of the static penetration length. 
This surface current exerts a force on the vortices, deforming the vortex array. Flux-line 
interactions cause the surface deformation to propagate in the interior of the sample, while 
impurity pinning and friction tend to impede its propagation. Using standard descriptions 
of flux-line dynamics, one finds that for a periodic perturbation modulated at frequency 
u, the penetration of the applied field is governed by a single complex penetration length 
A ac (o;) [1,8,9]. This length is closely related to the skin depth in normal metals and 
determines directly the surface impedance and the ac permeability of the material. The 
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peak in the ac dissipation occurs when X ac (u) is of the order of the sample size and is 
conventionally interpreted as associated with the thermally activated depinning of the flux 
lines More recently the peak in the ac dissipation in disordered superconductors 

has been interpreted as the signature of a continuous transition from a high temperature 



flux-line liquid to a low temperature vortex glass ||12|| . Good fits of ac measurements in 
single crystals and in thin films of YBCO have been obtained with the scaling relations 
of the phenomenlogical vortex glass model @J3|. In contrast the dissipation peak from 
mechanical oscillator data in clean samples has been interpreted in the context of a melting 
of the flux-line crystal to a flux-line liquid [^]|| . 

On the other hand, in a large part of the region of the (H, T) phase diagram probed 
by both dc and ac transport experiments the flux-line array is in a liquid state. In this case 
it clearly no longer makes sense to invoke the traditional description of flux-line dynamics 
|T^] in terms of the collective motion of crystalline flux bundles. The question then arises 
of how to incorporate collective effects in the description of flux-line dynamics in the liquid 
phase. 

A popular phenomenological approach for describing the dynamics of flux-line liquids 
is the " thermally activated flux-flow " ( TAFF ) idea [P, pTl|Jl^] , which essentially assumes 
an ideal gas of disconnected flux elements moving in a tilted washboard potential. This 
method generalizes the usual Bardeen-Stephen (BS) flux flow model [15[ to incorporate 
thermally- activated depinning of the flux lines. At temperatures higher than a typical 
pinning energy it reduces to the BS model. The TAFF phenomenology has been used to 



describe both dc |If| and ac response ||. This picture takes, however, little account of 
intervortex interactions in the flux liquid and is not very useful for describing a system 
undergoing a phase transition, where collective effects can lead to very large or diverging 
correlation lengths near the transition. 

As discussed elsewhere, a useful framework for the description of the long wavelength 
dynamical response of a flux-line liquid is hydrodynamics [|17| . In this context intervortex 
interactions are naturally incorporated in a viscosity that increases and goes to infinity 
at the phase transition from a vortex liquid to a vortex crystal or glass. In fact the be- 
havior of the viscosity at the transition serves as a signature of the order and nature of 
the transition itself. The hydrodynamic model naturally incorporates the intrinsic nonlo- 
cality of the electro dynamic response of a viscous flux liquid. Marchetti and Nelson used 
hydrodynamics to describe the response to dc perturbations [ |l8yr?|| . They showed that a 
large viscosity, as one can have as the phase transition is approached from the liquid side, 
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allows the effect of large scale spatial irihomogeneities, such as macroscopic twin bound- 
aries, to propagate large distances into the flux liquid and slows down considerably the 
motion of the flux lines. More recently Huse and Majumdar [|19| used the hydrodynamic 
model to describe multiterminal dc transport measurements in the vortex liquid regime of 
the cuprate superconductors. By calculating the dc response of a viscous flux liquid for 
realistic geometries including the effect of the sample boundaries, they demonstrated that 
the nonlocal resistivity arising from a finite flux liquid viscosity can explain the voltage 
patterns seen in experiments |2(| . 

In this paper we use the hydrodynamic model to describe the response of a viscous 
flux-line liquid to an ac perturbation applied at the surface of the sample. A large viscosity 
can impede the penetration of a surface field, since the viscous flux liquid cannot quickly 
adjust to the external field change. This is reflected in a nonlocal relationship between 
the vortex liquid flow velocity, which determines the electric field from flux flow, and 
the applied fields and currents. The nonlocality arises from viscous forces, which are 
proportional to the second spatial derivatives of the vortex flow velocity, and incorporates 
the force that remote fluid elements exerts on each other via intervortex interactions and 
entanglement. As a consequence of this nonlocality, the amplitude of the penetrating ac 
field is given by the linear superposition of two exponentially decaying contributions with 
different frequency-dependent penetration lengths, Ai and A2. Recently Sonin et al. [pT 



showed that the penetration an ac surface field that generates a periodic tilt of the flux 
lines is also governed by two characteristic length scales. In this case the new penetration 
length is associated with an elastic tilt mode of the vortex array that can impede surface 
penetration and was neglected in previous theoretical studies. In a viscous flux liquid for 
surface perturbations that compress or shear the vortex array there is a viscous mode 
that incorporates collective effects from intervortex interactions and entanglement and is 
responsible for the additional surface screening. This is to be contrasted to conventional 
TAFF or flux flow models, where the resistivity is local and the penetrating field has the 
form of a single exponential controlled by the single length A ac (w). The largest of the two 
penetration lengths, Ai(w), controls field penetration in the bulk of the material. At low 
frequency it simply describes Meissner screening. When the fluid viscosity is negligible 
\i(u) ~ X ac (uj). The additional penetration length, A2(w), is associated with viscous 
screening that prevents the magnetic field from building up to its maximum value right at 
the surface. The maximum of the penetrating field is then found at a distance of the order 
of |A2(w)| inside the sample. Viscous screening is effective in a thin surface layer of width 
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|A2(w)| ( I A2 (^) I << |Ai(w)| at all frequencies) and leads to the nonlocality of the electric 
field from flux motion. In a dense vortex array at low frequency | A2 (^) | is of the order of 
the intervortex spacing. 

In Section 2 we discuss the basic equations needed to describe the electro dynamic 
response of a type-II superconductor in the mixed state for the case when the flux array 
is in a liquid regime. These are the usual coupled Maxwell and London equations and 
a set of hydrodynamic equations for the viscous flux liquid. In Section 3 we present the 
solution for an ac field penetrating at the surface of a semiinfinite superconductor. The 
role of viscosity and its effect on macroscopic properties, such as the material's surface 
impedance are discussed. The solution for a superconducting slab is given in Appendix A. 



2. Flux-Line Hydrodynamics 

We consider a uniaxial type-II superconductor in a static field H applied along the 
c axis of the material, which is chosen as the z direction, i.e., H = Hz. The static field 
produces an array of flux lines that are on the average aligned with the field direction. The 
mean spacing ao of the flux lines in the xy plane is determined by the average equilibrium 
induction B = B z that penetrates the sample, with ao — l/y/n Q and n = -Bo/Vo 
the areal density of flux lines (<j)Q = hc/2e is the flux quantum). For H » Hqi and 
neglecting demagnetizing effects, Bq ~ H. For high k superconductors, such as the HTC 
materials, there is a large part of the (if, T) phase diagram where << ao << \ a b (i.e., 
Hci « B « Hci), with A a b and ^ a b the penetration and coherence lengths in the ab 
plane, respectively. Here we restrict ourselves to this region. In addition, we are interested 
in properties of the flux array on length scales larger than the mean intervortex spacing 
where one can describe the system in terms of a few hydrodynamic fields coarse-grained 
over several flux-line spacings. 

The electrodynamics of a type-II superconductor is described by Maxwell's equations 
for the local fields b(r,t) and e(r, t), 

Vxe + V = 0, (2.1) 

c 

„ , 4tt . 1 „ , 

V x b = — j + -d t e, (2.2) 

c c 

with V • e = 47rp and V ■ b = f2"2]j . These equations have to be supplemented with 
a constitutive relation for the current density, j. Following Coffey and Clem H], we use 
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a two fluid model and write j = j n + j s . The normal current density j n is specified 
by Ohm's law, j n = a n e, with a n the electrical conductivity of the normal metal. The 
constitutive equation for the supercurrent j s is obtained by minimizing the Ginzburg- 
Landau free energy functional for an anisotropic superconductor in the presence of vortices. 
As customary, the anisotropy is incorporated in an effective mass tensor. In the London 
approximation, where the magnitude of the order parameter is assumed constant and only 
fluctuations in the phase 9 of the superconducting order parameter are retained, we obtain 

A.j s = -^-(A-0 o V0), (2.3) 

where A is the vector potential, with b = V x A, and A is a diagonal tensor with com- 
ponents A xx = A yy = \ 2 ah and A zz = \ 2 C . Here A c is the penetration length along the c 
direction, with X c /X a b = ^/m c /m a b = p » 1. Taking the curl of Eq. ( |2.3|) and averaging 
over lengths large compared to the intervortex spacing, one obtains the London equation 
in the presence of vortices, 

VxA-j s = £-(-b + ^ T). (2.4) 

Here T = zn + r, with n(r, t) and r(r,t) the coarse-grained hydrodynamic density and tilt 
fields of the flux array, respectively. The microscopic density and tilt fields are given by 



/) 

i=l 



{r,t) = Y;6{r ± -r i {z,t)), (2.5) 



and 

N 



r mic (r, = £ ^f^S(r± - r,(z, t)), (2.6) 

i=l 

where Ti(z) is the position of the z-th vortex line in the xy plane as it wanders along the z 
direction, and r = (rj_, z). The tilt field describes the local deviation of an element of flux- 
line liquid from the alignment with the z direction. The hydrodynamic fields are obtained 
by coarse graining the corresponding microscopic fields over several vortex spacings. The 
density and tilt fields are not independent, but are related by the constraint that flux-lines 
cannot start nor stop inside the medium. This requires 

d z n + V±-r = 0. (2.7) 
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Using Maxwell's equations one can eliminate the current from Eq. ( j2.4|) in favor of 
the fields b and e. This can be done by multiplying Eq. (|2.2| ) with the tensor A from the 
left and then taking the curl, with the result, 
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VxA-(Vxb)+b = (P T-^(X 2 ab d t b-X 2 c V ± xze z )--d t (^d t b-X 2 c V ± xze z ). 



The second term on the right hand side of Eq. ( |2.8| ) comes from the normal part of the 
current density. The last term arises from the displacement current and is negligible for 
all frequency uj « % ~ 10 15 Hz. Since we are interested in experiments carried out at 
frequencies no higher than microwave, we will neglect this term here, even though it can 
be easily incorporated in the calculation. 

To describe the response of the superconductor to a perturbing field 5H a applied 
at the surface of the sample, one needs to solve Eq. (|2.8|) with appropriate boundary 
conditions on the field. The local field B is, however, coupled to the vortex distribution 
through the right hand side of Eq. (|2.8| ) and the vortices in turn move in response to 
changes in the field, as discussed for instance by Coffey and Clem ||. Equation ( |2.8|) has 
to be supplemented with equations describing the vortex dynamics. Here is where our work 
differs from that of other authors. We assume the flux-line array is in a liquid state and 
describe its dynamics through the set of hydrodynamic equations for the density and tilt 
field discussed elsewhere. The hydrodynamic equations consists of a continuity equation 
for the density, 

d t n + V±-j v = 0, (2.9) 
with j„ = nv, and a continuity equation for the tangent field, 

d t r a + dpj T aP = d z j va , (2.10) 

where j^g are the components of the tangent flux tensor. This is a 2 x 2 antisymmetric 
tensor and has therefore only one independent component, j T xy = € a gj^ g = —jy X - At equi- 
librium n = no and the tangent field vanishes, To = 0, since the flux lines are on the average 
aligned with the applied field. In the following we will only discuss the linear response 
of the flux liquid to external perturbations of small amplitude [23[]. The hydrodynamic 



equations can then be linearized in the deviations of the hydrodynamic fields from their 
equilibrium values, Sn = n — no and St — r, with j„ rs nov, 

d t 8n + n Vi • v = 0, (2.11) 
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d t r a + dpjlg = n Q d z v a . (2.12) 



As discussed in [|17| , these equations have to be closed with constitutive equations for the 
fluxes. Neglecting the Hall current, and keeping only terms linear in the fluctuations from 
equilibrium, the constitutive equations for j,, and are given by, 

-TV + ?? s Viv + //, ( Vx(V ± ■ v) + n z 8 2 z w - -B x j = 0, (2.13) 

c 

-lrj T xy + r} T Vlf xy + Vr Z d 2 j: y + ^B • j = 0, (2.14) 

where j = (c/4-7r)V x <5B, with SB = B — Bo- The first terms in both Eqs. ( |2.13j ) and 
( [2.14] ) describe the drag on the flux lines arising from interaction with the crystal lattice, 
with t arid Tr friction coefficients per unit volume. The contribution to t from interaction 
of the normal core electrons with the underlying crystal lattice can be approximated by 
the Bardeen-Stephen coefficient ^bs{T,H) = no7rfoV n /2e 2 ^ b |15|]. Weak point pinning 



centers can be approximately incorporated in a renormalized friction by assuming t — 
r )BS eUp ^ kBT i where U p is a typical pinning energy, as done in thermally activated flux 



flow (TAFF) models ||14|1 . The friction coefficient Tr can be related to the relaxation rate 
of an overdamped helicon and we estimate Tt ~ T- The next three terms in Eq. ( |2.13j ) 
and the next two terms in Eq. (|2.14|) describe the viscous drag in the flux-line liquid 



from intervortex interactions and entanglement. These effects are incorporated in the 
viscosity coefficients: rj s , n^ and r] T , denoting the shear and bulk viscosity coefficients of an 
anisotropic liquid, and the "tilt" viscosities n z and rj TZ associated with velocity gradients 
in the direction of the applied field. Finally, the last two terms on the left hand side of 
both Eq. (|2.13| ) and Eq. ( |2.14|) represent reversible forces on an element of flux-line liquid. 



To linear order in the hydrodynamic densities these can be written as 

--B xj = -— / dr'[Vic L (r-r>n(r')-a z cJ 4 (r-r')r(r')], (2.15) 
c n J 

^B -j = J dr'c^(r-r')Z-[V ± ' xr(r')]. (2.16) 

Here cl(t) and c^(r) are real space elastic constants. Their Fourier transforms, c^(q) and 
€44 (q), are the nonlocal compressional and tilt moduli, respectively |Tj]]. The liquid tilt 
modulus can differ from that of an elastic flux array because of flux-line cutting. We note 
that the nonlocality of the elastic constants is automatically incorporated in our treatment. 
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When discussing the response of a viscous vortex liquid to an ac perturbation one 
needs to modify Eqs. ( pT9| ) -( |2.14j) to include viscoelastic effects analogous to those present 



in entangled polymer melts ||24j| . These effects will be important in the vortex liquid in 
the region of temperatures just above the transition to a solid phase because the viscosity 
gets very large in this region. If velocity gradients are small, viscoelastic effects can be 
incorporated following the theory of polymer dynamics by modifying Eqs. ( |2.13| ) and ( |2.14| ) 

as 

- 1 v ce (r,t) + dt' J cir / C a ^(r-r / ,t-t / )^^(r / ,t / )-^[B xj(r,t)] a = 0, (2.17) 
and 

-7rj^(r,t) + J*dtfJ dv'CJjir -v',t- t^d'^v' ,t') + ^B • j(r,t) = 0, (2.18) 

where latin indices run over the three cartesian coordinates, % = x,y, z, while Greek indices 
only run over the two coordinates normal to the applied field, a = x,y. The components 
of the tensors C a ipj(r,t) and CjAr,t) are the Green-Kubo correlation functions that de- 
termine the viscosity coefficients. It is convenient to introduce the Fourier and Laplace 
transforms of these correlation functions, defined as 

C ai()j {q,u) = J dt J dre ifi - r C aWj (r,t) (2.19) 

for Im{uj) < 0. There are five independent Green-Kubo correlation functions, defined by 

C ai f3j(q,u}) =f) s {(\,u)5 a p[5ij -S iz S jz ] + fj z (q_, io)5 a p5 iz 5 jz 

1 (2-20) 
+ - [r)i(q,u) - fj 8 (ci,uj)][8 a iSi3j + SajSpi], 

and 

Clj(q,u) = fj T (q,u)8 a /3[Sij - S iz S jz ] + fj TZ {(\,u})5 a p5 iz 5j z , (2.21) 

where ?)(q, u) for v = s,z,1,t,tz are the frequency- and wavevector-dependent viscosity 
coefficients. These determine the static viscosity coefficients according to, 

Vu =?)(q = 0,w = 0), (2.22) 
for v = s,z, /, r, rz, where rji = r] s + Vb is the longitudinal viscosity. 
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In entangled polymers the relaxation of shear and compressional stresses described by 
the Green-Kubo correlation functions is nonexponential |[25|| . This reflects the existence 
of a wide distribution of relaxation times in the systems. The same feature is expected to 
occur in strongly interacting flux-line liquids. Here for simplicity we model the time decay 
of the Green-Kubo correlation functions as exponential according to the phenomenological 
Maxwell model of viscoelasticity [^6|. The wave vector and frequency-dependent viscosities 
are then given by 

VM^)= ^ (2.23) 
1 + uvT u (q) 

The characteristic relaxation times r v are chosen so that the viscoelastic model incorporates 
the essential feature that the fluid behaves like a viscous liquid on long times scales and 
as an elastic solid on short time scales. In other words for ui » 1/t v for all z^, the 
hydrodynamic equations should reduce to the equations of continuum elasticity for a flux- 
line array. This requires 

r s = Vs/G, (2.24) 
7* = W(cii(q) + G-c L (q)), (2.25) 
T 2 =W( C4 4(q)-c;L(q)), (2-26) 

where cu, G and C44 are the compressional, shear and tilt elastic moduli of a flux-line elastic 
medium. For an Abrikosov flux-line lattice G ~ cqq and the wavevector dependence of the 
shear modulus cqq is always negligible. The compressional modulus is given by cn(q) = 
cj,(q) + C66 and Cn(q) ~ cj,(q) since cqq « cj,(q) at all but the largest wavevector 



q ~ ksz with ksz = a/ 4-7r/ao the Brillouin Zone boundary, where cl^bz) ~ cqq. As a 
consequence we can approximate cu by cl in eq. (2.25) with the result r; ~ rji/G. The 
wavevector dependence of both the shear and compressional relaxation time is therefore 
negligible. 



3. Surface Impedance and AC Permeability 

In this section we discuss the response of a flux-line liquid to a weak ac field 5H a = 
iSHae 1 ^ 1 applied parallel to the static field H, with SH a « H. For clarity we first 
consider the case of a semiinfinite superconductor occupying the half space y > 0. The 
perturbing field is applied at the surface, y = 0, and generates a surface current in the 
+x direction. The surface current in turn exerts a Lorentz force normal to the sample 
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boundary that yields a compression of the flux array. To evaluate the change in induction 
in the superconductor as a result of this perturbation we need to solve the coupled set of 
Maxwell and London equations and the hydrodynamic equations of the flux liquid. Since 
the lines remain on the average aligned with the z direction, the change in induction will be 
of the form 5B = z5B z (y)e' lujt and r = 0. The linearized hydrodynamic equations reduce 
to 

dv 

iuj5n(y) + n —^- = 0, (3.1) 
dy 

-ivyiv) + mH^t - ^rUv) = o, (3.2) 

where j x (y) = (c/47r)(dSB z /dy) is the total current density, including the response of the 
medium. The field 5B includes both the Meissner response to the applied field and the 
change in induction generated by the motion of the vortices. It is determined by the 
London equation (|2.8|) . For the simple geometry considered here this becomes, 

=&*•<»>. <»■»> 

where A n / = a/ c 2 /4niua n is the normal fluid skin depth and A is the London penetration 
length in the ah plane. For simplicity we drop in this section the subscript on A a &, since 
this is the only penetration length that is relevant for the geometry considered. Equation 
( p73|) has to be solved with the boundary condition 5B z (y = 0) = 5H a . 

When fji = 0, corresponding to a simple flux flow model in the absence of pinning, one 



can eliminate 5n between Eqs. ( |3.3| ), ( |3.1| ) and (|3.2|) to obtain a second order differential 
equation for 5B Z , which has solution 

5B z (y)=5H a e-y/ x ^\ (3.4) 

where 



AacM 



X 2 + X 2 f (u) n l/2 



(3.5) 



with \f(cu) = a/cl(0)/zw7 the ac penetration length of a flux liquid ^7|. As in a normal 



metal, the response of the material to an ac field is entirely characterized by a single length 
A ac (a>), that plays the role of the skin penetration depth in metals. The inclusion of the 
normal current guarantees that X ac (uj) reduces to the normal metal skin depth at H C 2 where 
A diverges . The frequency dependence of X ac is controlled by two frequency scales. The 
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first, ujf = c^(0)/7A 2 is associated with flux diffusion: 1 represents the time it takes for 
the field to diffuse a length A into the sample. The second frequency scale is the frequency 
governing ac field screening by the normal fluid component, u n f = c 2 /(4-Ka n X 2 ). If the 
friction coefficient 7 is approximated by its Bardeen- Stephen value, we find Uf/u n f = 
B / H C 2 < 1. The ratio of these two frequencies is even smaller if 7 is renormalized to 
incorporate the pinning by weak point defects as in TAFF models, since in that case 
u)f/u n f ~ (B/Hdz) exp(— Up/hsT). The frequency dependence of A ac (u;) is characterized 
by three regions, as shown in Fig. 1. For uj « Uf, X ac (uj) ~ A/(u;) ~ uj~ 1 / 2 and the 
field penetration is controlled by flux flow. For ojf«u« uj n f, A ac (u;) ~ A and one has 
essentially static Meissner screening since the flux lines cannot flow on experimental time 
scales. Finally, for to » u n f, \ ac (uj) ~ X n f(uj) ~ u~ x l 2 and the electrodynamic response 
is controlled by the normal fluid. We remark that the frequency dependence of \ ac {uj) in 
the region where penetration is controlled by flux flow is the same as in a normal metal. 
In the flux liquid model considered here, weak point disorder decreases Uf and widens the 
region of Meissner screening. As pointed out by Geshkenbeim et al [31]., iOf is very large 
(typically ujj ~ 10 14 Hz) and in the frequency range where the experiments are usually 
carried out (1 — 10 6 Hz) one has X ac ~ A/. 

The quantity that is usually measured at these frequencies is the surface impedance 
Z s (u) = R s + iX s , defined as 

E x (0) 4tt E x (0) 

Z s {u) = r oo . . r-, = — , nV (3.6) 
Jo 3x(y)dy c bB z {0) 

where ^(0) is the electric field at the surface. The surface resistance, R s , and the surface 
reactance, X s , determine the loss and the resonance frequency, respectively, when the 
superconductor forms part of a resonant circuit. When fji = the surface impedance is 
given by the familiar expression, 

Z s (u) = ^-^-\ ac (u). (3.7) 
c z 

It is related to the ac resistivity, p ac {uj) = E x (y)/j x (y), by Z s (u) = [4niujp ac (uj)/c 2 ] 1/2 . 
If I A/ 1 >> A, the ac resistivity is independent of frequency and it is simply given by the 
Bardeen- Stephen dc resistivity, p ac = Pf = 47tcl(0)/c 2 7 P7] . 

Finally, for a superconducting slab of thickness 2W occupying the region of space 
\y\ 5; W, one can use this simple flux-flow model that neglects intervortex interaction to 
evaluate the ac permeability, n(u), defined as 
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with the result, 

^-w^O- (3 - 9) 

The ac permeability depends on the sample thickness W through 

Aq C _ r \ 2 /W 2 +UJ D /iU ,l/2 

W 1 l + iu/u nf 1 [ ' ' 

where ujo = ujf\ 2 /W 2 is the inverse of the time for flux diffusion across the sample. The 
normal fluid contribution is negligible provided u « u n f, where u n f >> too- In this 
range of frequencies a crossover in fx and a peak in fx occur when |A oc | ~ W or u> ~ ud, 
i.e., the time modulation of the applied field matches the time for flux diffusion across the 
sample thickness. If u « u>d, \X ac \ » W and the field penetrates completely, yielding 
fx =1. If u » u>d the flux liquid cannot diffuse on the time scale over which the 
perturbation changes and the applied field can only penetrate a surface layer of width A, 
yielding fx « X/W tanh(W/ X) << 1. The crossover between these two limits is marked 
by a maximum in the absorption fx . Finally, for u ~ u n f the normal fluid contribution 
becomes important and leads to a second crossover from fx « X/W tanh(W/ X) for cud << 
uj « u> n f to fx ~ for oo » uj n f, accompanied by a second peak in fx . When u » u> n f 
the normal quasiparticles cannot follow the changes in the external field which is completely 
expelled from the sample. If W » A, both fx and fx are very small at u ~ w n /, where 
the normal fluid contribution becomes dominant. As a result, the second peak in fx is 
very small, as show in Fig. 2a, where W/X = 10. On the other hand, if W > X the second 
peak in fx at u ~ u n f can exceed the peak from flux flow ( Fig. 2b). In general, the 
frequencies u n f is, however, outside the range of frequency usually probed in experiments. 

Recently Coffey and Clem incorporated flux creep and pinning in the model described 
above by introducing a frequency-dependent single-vortex mobility (the inverse of our 
friction coefficient) ||. They found that, while pinning and creep affect the location of 
the peak in the absorption, no new length scale appears in the electrodynamic response, 
which is qualitatively similar to that obtained from simple flux-flow models. 

In a viscous flux liquid intervortex interactions and entanglement as described by the 
viscous force in Eq. ( |3.2| ) make the electrodynamic response nonlocal. This nonlocality is 
governed by the new viscous length 



6(u) = 

12 



(3.11) 



that controls the relative importance of the first two terms of the left hand side of Eq. 
(3.2). To clarify the meaning of this new length scale it is useful to eliminate the flux 
density and flow velocity from Eq. (3.1)-(3.3) in favor of the electrodynamic fields 5B z (y) 
and E x (y) = ^v x (y) + ^^j s , with the result, 

x ~ ~W = pf3xiy) + V' ~ "V ^ ' ( } 

where Pf = Bq/c 2 ^ is the flux flow resistivity, and 

XV ^-(,L + ^+^ = 0. (3,3) 



The second term on the right hand side of ( |3.12| ) arises from Meissner screening of the 
fields over a surface layer of width A. To understand the physical meaning of the new 
length scale 5(u) we define a contribution from flux motion , E[, to the electric field by 
E x = E[ + AltX l iu) j s . Substituting this into , we obtain 

El-8 2 ^f=p f jM. (3.14) 



When the viscosity is very small the second term on the left hand side of Eq. (|3.14|) is 
always negligible and the field from flux motion is simply determined by Ohm's law for 
a normal metal of resistivity given by the flux flow resistivity, pf. The resistive response 
of the medium is local and the ac resistivity p ac is independent of the frequency, with 
Pac = Pf- When 6 is sufficiently large the resistive response of the medium is nonlocal 
: the electric field at r is determined by the current at spatially remote points r' as a 
result of the force that remote fluid elements can exert on each other via interactions and 
entanglement. 

Using the Maxwell model of viscoelasticity given in Eq. (2.23), the frequency depen- 



dent viscous length 5(uj) is given by 5(uj) = Sq(1 + iuri) -1 ! 2 where Sq = \fr\ip) is the 
static viscous length discussed earlier by Marchetti and Nelson fl7[j . The static viscosity 
of a flux-line liquid has been estimated elsewhere employing analogies with the physics of 
entangled polymer melts |T7| , |28[] . Assuming 7 w 7^5, one finds do ~ ao exp(C/ x /2/csT), 
where U x (H, T) is the typical energy barrier for flux-line cutting, which is expected to 
vanish at H C 2- If the barriers to flux cutting are sufficiently large, at low temperatures 
the vortex array can get stuck in a polymer glass regime of entangled lines, character- 
ized by infinite viscosity on experimental time scales. A simple estimate of the crossing 
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energy gives U x ~ 2(y/2 — 1) , /^^-aoeo In /?, with eo = (<^o/47rA a f,) 2 p9| . This can also 



,/ 2 c„«„ / m ab 



be written in terms of the clean flux lattice melting temperature T m = aL L e$ao, 
with aL ~ 0.15 ~ 0.3 the Lindemann parameter J29|, as U x /T m = c x /a|, with 
c x = 2(y / 2 — l)lnK. This shows that the crossing barrier associated with entanglement 
can become very large above T m and preclude crystallization on experimental time scales. 



Recent calculations of the energy U x by Moore and Wilkin |3(| and by Carraro and Fisher 
1 31] confirm these simple estimates. Using ctL ~ 0.3 and k ~ 200, we find 5q ~ aoe 50Tm ^ T . 
It is then clear that the static viscous length 5q can become very large in the flux liquid 
regime and the resulting nonlocality of the dc response can be probed experimentally. 

Viscous effects introduce two new frequency scales in the ac response of a flux array. 
The first is the frequency ui = 1/r/ = G/rji describing the relaxation of shear stresses and 
controlling the frequency dependence of the viscous length 5(uj). The second frequency 
scale is defined by |A/(o>)| ~ 5q, corresponding to u ~ u v = cl(0) /r)i. Both ivi and u v 
decrease with increasing static viscosity and in flux arrays u>i « u v since G << c^(0), 
as discussed below. The frequency u)\ only enters through the frequency dependence of 
S(u) and governs the crossover from liquid-like to solid-like response of the flux array as a 
function of the frequency uj of the external probe. To understand the crossover it is useful 
to first neglect the fourth order derivate of the field in ( 3.13 ) . The single length scale 
governing ac field penetration is then A^ c + (5 2 = X 2 + Xj + S 2 . At low frequency (u << ui), 
5(u) w Sq and X 2 C w A 2 + X 2 + 5q, with A/ = a/cl(0)7^7- The model then simply 
describes penetration via flux diffusion (flux flow) of a vortex liquid of static viscosity r\i. 
In this case the viscous force provides an additional static damping of the penetration field. 
This additional damping is, however, negligible if u < w v . If G « cc(0), then u>i « uj v 
and the flux flow contribution to the penetration length always dominates the contribution 
from 5q in this low frequency regime. At high frequency ( u » ui ), 8 (u) ~ a/ G/iwy and 
X 2 + S 2 ~ A 2 + e f W _)_ -Q_ _ y2 _|_ c n(°) m j n ^^jg regime field penetration is governed by 
flux flow of a vortex lattice. The model therefore describes the crossover from liquid-like 
to solid-like response at the frequency This frequency decreases as the static viscosity 
increases and uii/iVf = (G/cl)(X 2 /Sq). Assuming G ~ cqq, where cqq is the shear modulus 
of the Abrikosov flux lattice, we find G/c L (0) w (1 - 6) 2 /4k 2 6 2 , if b = B/H c2 > 0.25, and 
c 6 6/cl(0) w 1/4k 2 6 2 if 0.3/k 2 < b < 0.25. ^ Using material parameters of YBCO at 



85K, we estimate G/cl(0) < 10 ||32] , p3l . As a result, at the frequency uj ~ u>i ( 



when the crossover from liquid-like to solid-like behavior occurs ) the conventional flux flow 
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contribution to the penetration length is so large that it always dominates the ac response 
of the flux array. This result is not unexpected. It simply reflects the fact that the only 
difference in response between a vortex solid and a vortex liquid arises from the small 
difference in the compressional moduli. A perturbation that generates a compression of 
the flux array is therefore not a good probe to distinguish between liquid-like and solid-like 
response. Viscous effects may be observable only provided G/cl ~ 1, as discussed below. 

We now discuss in detail the field penetration for the semiinfinite geometry when 
fji 7^ 0. It is clear from (|3.13|) that an additional boundary condition is needed to solve 
the equations. The additional boundary condition used here is Sn(y = 0) = 0, which 
follows from the analysis of the fields generated by the vortices and their images near the 
surface |T0[1 . Assuming that all perturbations vanish as y — > oo, we obtain 



8n(y) 



\2 



A 



A| 



-y/M 



(3.15) 



where X±(uj) and X2(uj) are two frequency-dependent complex penetration lengths, given 
by, 



A?» 
The fields are given by 

sbm 5Ha 

and 
EM = 1 



A) 



A| 



1 

A^ 



4X 2 5 2 



(al + 52)2 ~TTava 2 



nf 



(3.16) 



\ 2 



-V/M 



1 



1 



\ 2 



(3.17) 



■iuj 



5H n 



( \2 \2 ^ 



A| 



A 5 



\ 2 



)A ie -^-(^-^ + ^)A 2 e-^ 



(3.18) 



In a viscous flux liquid the penetration of the ac field is governed by two length scales, 
Xi(iv) and X^ico)- When fji = 0, A2(w) = and Xi(uj) = X ac (uj). In this case Eq. ( |3.17p 
simply reduces to (|3.4|) . The role of the new penetration length, A2, can be understood by 



examining the change 5n(y) in the vortex density arising from the ac field, given in ( 3.15 



The density 8n(y) reaches its maximum near y Q w [| Ai 1 1 A2 1 / ( | Ai 
I Ai| >> I A2 1 for all frequencies of interest here, yo 
The density grows from zero at y = to its maximum value at yo 
to zero over a length | A 1 1 >> | A2 1 - I n other words both the surface currents associated 
with the Meissner response (A 7^ 0) and the spatially inhomogeneities in the electric field 



|A 2 |)]]n(|Ai|/|A 2 |). Since 
|A 2 | and Sn(y ) « (5H a /(f) )X 2 f /X 2 . 

I A2 1 and then decays 
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arising from intervortex interaction (fji 7^ 0) impede the build up of an appreciable vortex 
density in a surface layer of width | A.2 1 . When either fji = or A = the width of this 
surface layer vanishes. A model that neglects variations on the length scale of order A 
(this corresponds to letting A = in Eqs. (3.13-16)) assumes that 4>odn(y) = 5B z (y) 
and neglects the surface currents responsible for the jump in the tangential component 
of the average magnetic field at the surface. Similarly, when 5 = one neglects spatial 
inhomogeneities in the electric field from flux motion at the surface and the flux-line density 
again has an unphysical finite value at the surface, given by 5n(0) = (SH a /(f) )X'j/(X 2 j: + X 2 ). 
It is also possible to satisify the boundary condition 5n(0) = in a microscopic model 
that incorporates the interactions of the vortices among themselves and with their images. 
This boundary condition cannot, however, be satisified in the conventional flux flow model. 
The hydrodynamic model described here provides a phenomenlogical description of flux 
dynamics that can satisfy this boundary condition because it incorporates the nonlocalities 
in the velocity and electric fields due to intervortex interaction. 

Due to the nonlocality induced by the viscous screening the electrodynamics response 
of the system needs to be described in terms of nonlocal response functions. The local real 
space ac resistivity and ac permeability are defined by 

*to) = /<V/U» -✓,«».(✓). (mo) 

and 

5B z (y) = J dy'^ ac (y - y', uj)H a {y'). (3.20) 

In an infinite system the above nonlocal relationships in real space simply yield the usual 
linear relationships between the Fourier components of the fields and currents, 

E x (q, u) = p ac (q, w)j x (q, u) (3.21) 

5B z (q, u) = jU oc (q, v)H a (q, u) (3.22) 

Due to this nonlocality the ac resistivity is not simply related to the surface impedance. 
The voltage drop measured in a transport experiment has to be evaluated for each specific 
experimental geometry using realistic boundary condition, as discussed in ||19[| . The effect 



of the viscosity and the crossover between liquid-like and solid-like behavior should be 



observable in ac multiterminal experiments of the type described in [ 2"0[] . Similarily, the 



nonlocality of the permeability could be probed by measuring local magnetization profiles 
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inside the sample. The net field penetration in a slab of thickness 2W is still described by 
the macroscopic permeability defined in ( |3.8|) , which corresponds to the q = component 
of //(q, u). Its expression is given in Eq. (A. 2). 

To simplify the discussion of the two length scales Ai(o;) and A2(w) and of other 
electro dynamic properties we will drop the normal fluid contribution below. As in the case 
r\i = 0, the normal fluid contribution introduces an extra crossover at very high frequency 
(u ~ u n f ~ 10 15 Hz) and can be easily distinguished from all the other relaxational modes 
which occur at frequency scale at least one order of magnitude less than u n f. 

The new length scale A2 is in magnitude at most of the order of A at all frequencies. 
When the static viscosity is small, i.e. So/X << 1, A2 is always negligible compared to 
both Ai and A, as shown in Figs. 3a and 4a. In this case the field penetration is governed 
by the single length scale Ai, with Ai ~ A/ over the entire frequency range of interest ||34j| . 
The ac permeability is well approximated by the familiar formula (jl(uj) ~ Ai(u;)/Wtaiih 
(W/Xi(u)) and is shown in Fig. 3b. Experiments are typically carried out at frequencies 
u < u>f. As discussed earlier, the transition from solid-like to liquid like response takes 
place at u ~ Since u>i « Uf in a viscous liquid due to the small value of the ratio 
G/cl, the crossover at to ~ uji occurs well into the flux flow regime, in the sense that 
I Xf (lu) I » Sq and conventional flux flow dominates the response, as shown in Fig. 3b. 

The expression used above for the shear relaxation time 77 is simply a phenomenlogical 
estimate. In particular the moduli G and cl(0) appropriate for an entangled polymer 
glass are not known. An alternative approach would be to simply consider t\ or u)\ as a 
parameter. The effect of viscosity on the ac permeability will then be observable if the 
G/cl is not too small, while the viscosity is sufficient large. To illustrate this we show 
in Fig. 4b-6b the ac permeability for ui/uf = 0.01 and Sq/X = 1,10,100. When the 
viscosity is small (Sq/X = 1), A2 is still negligible at all frequencies and the response is 
given by conventional flux flow with the maximum at \x corresponding to |A/| ~ W (Fig. 
4b). When Sq/X » 1, A2 can become of order of A and the width of the surface layer 
where flux penetration is impeded is no longer negligible. In this regime and for uj < uj\ 
we can approximate X\ ~ X\ c + 5q and A2 ~ 2A 2 5q/(A^ c + Sq). There is a new crossover 
at the characteristic frequency uj^ where |A/| ~ 5q , as shown in Figs. 5a and 6a. For 
uj « uj v we find Ai ~ A/ while A2 ~ Sq^wJuj} is negligible. For uj » uj v Ai ~ 5o 
and A2 ~ A. Correspondingly, a second peak develops in /i at u ~ for intermediate 
viscosity (Fig. 5b). At very large viscosity, the crossover in the ac response occurs at 
uj v . To understand this we recall that | A2 1 is the distance from the surface over which the 
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vortex density builds up to its maximum value, while |Ai| characterizes the decay of the 
density from its maximum value to zero within the sample. If uj v « uj\ ( as it is the case 
for the parameters of Figs. 6 ), |Ai| >> 2W at all frequencies (Ai ~ A/ ~ 1/-\/uj for u> < uj v 
and Ai ~ do for < u < u{) and the ac perturbation simply cannot penetrate into the 
sample, unless | A2 1 is essentially zero, as it is the case for u « u^. The real part of 
the ac permeability drops therefore from the value corresponding to complete penetration 
(fx = 1) to Meissner response (/x = 0.1) at u = u> v , where A2 becomes comparable to A 
and the viscous screening discussed above becomes appreciable. 

Finally, the surface impedance is obtained by inserting Q3.18|) in ( j3DD ,with the result 

C 1 (A1 + A2) 

The surface impedance for a superconducting slab of finite thickness 2W in the y direction 
can be calculated in a similar way. The result is given in Appendix A. Again, when 
IA2I << A << I Ai| the field penetration is governed by the longest length scale and Eq. 
(A. 3) can be approximated by 



ZJu) 



3-, 



Aitanh(W7Ai) + ^ 



^ftitjj , , r — — — 1 - , 

-Aitanh(W/Ai) 



(3.24) 



c 2 



4. Summary 

We have studied the linear response of a viscous flux-line liquid to ac perturbations 
by using a hydro dynamic theory. The flux array is described as a viscoelastic medium that 
responds elastically to perturbations varying on time scales shorter than the characteristic 
time u^ 1 for relaxation of shear stresses, while it behaves as a viscous fluid on time scales 
large compare to uif 1 . For realistic values of parameters for YBCO, the characteristic 
frequency uj\ is, however, small compared to the flux flow frequency Uf. As a result, 
the crossover from liquid-like to solid-like behavior occurs well into the flux flow regime 
and it generally cannot be detected in the ac permeability which probes the spatially 
averaged response of the system. This is because an ac permeability measurement probes 
the response of the flux array to a compression. Both flux liquid and flux solid have 
nonvanishing compressional moduli and the values of these moduli in the two regimes are 
very similar, in virtue of the very small shear modulus of the vortex lattice. A signature of 
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the viscous flux-line liquid that distinguishes it from a flux lattice is the intrinsic nonlocality 
of the response. This is responsible for the appearance of the second penetration length, 
A2. It should be possible to probe this nonlocality by flux imaging experiments, provided 
the spatial resolution is sufficiently high to detect variations of length scales of order of 
IA2I ~ A. 

This work was supported by the National Science Foundation through grants No. 
DMR-9112330 and DMR-9217284. 



Appendix A. AC response for slab geometry 



In this Appendix we describe the response of a superconducting slab occupying the 
region \y\ < W to an ac field zSH a e' lujt applied at the sample surfaces. The local magnetic 
induction in the superconductor is given by, 



SB z (y) 



5H n 



te- 



ll 1 . cosh(y/Ai) .1 1 1 cosh(y/X 2 ) 



■A| A 2 A2/cosh(W7Ai) y \j A 2 X 2 n /cosh(W/X 2 



(A.1) 



The ac permeability defined in (|3.8| ) is 



mM = (f) (T^[(^-^-4) Aitanh(w/Ai, "(^"^"i) A2tanhWA2) 

(A.2) 



The surface impedance at one of the boundaries is given by 
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~ A 2 
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~ A 2 
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A 1 tanh(W/A 1 ; 
A 2 tanh(W/A 2 ; 



(A.3) 
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Figure Captions 

Fig. 1 

Modulus of the complex ac penetration length as a function of reduced frequency uj/ujf for 
the zero viscosity flux liquid with uj n f /ujf = 50. All length scales are measured in units of 
the London penetration length. Three regimes can be identified (i) conventional flux flow 
(u << LOf, X ac ~ A/), (ii) static Meissner screening (u>f << u « w n f, X ac ~ A) and (iii) 
normal fluid dominated regime ( u » w n f, ^ac ~ A n /). 

Fig. 2 

Real and imaginary parts of the ac permeability for r\ = and u> n f/u>f = 50: (a) W/X = 
10.0 and (b) W/X = 2. 

Fig. 3 

AC penetration lengths (a) and permeability (b) as functions of the reduced frequency uj/ujf 
for a viscous flux liquid, with G/cl = 10~ 3 . No appreciable change in the ac response is 
observable for S 2 /X 2 = 1, 10, 100. The normal fluid contribution has been neglected here 
and in Fig. 4, 5 and 6 below. As a consequence fx — > ^ tanh(VF/A) as u — > oo. 

Fig. 4 

AC penetration lengths (a) and real and imaginary parts of the complex ac permeability 
(b) as functions of the reduced frequency uj/ujf for u v /ujf = 0.01 and Sq/X 2 = 1. The peak 
for /i at u! ~ u!f = 1 corresponds to |A/| ~ W. No effect of the viscosity is observable for 
S /X = 1. 

Fig. 5 

Same as Fig. 4 for S /X = 10. The incipient second peak which occurs at uj ~ uj v ~ 
cu/(A 2 /5o) represents the transition from viscous to conventional flux flow. This second 
peak is not present in Fig. 4b because the viscosity is too low there and the two charac- 
teristic frequencies and Uf coincide. 

Fig. 6 

s 2 

Same as Fig. 4 for So/X = 100. The crossover occurs at uj = ujq = u)f /(-$■) corresponds to 
the impediment of the penetration of ac perturbative field due to the viscous screening. 
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